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I. INTRODUCTION 



Recently there have been great interests in generahzed parton distributions (GPDs) of 
hadrons (see e.g. Refs. [1-3] for the original work, and Refs. [4-6] for recent reviews). GPDs 
relate quite different physical quantities, such as Feynman's parton distribution functions 
(PDFs) and hadron form factors, in the same framework. By generalizing PDFs' one di- 
mensional parton distribution pictures, GPDs provide the three-dimensional pictures [7, 9]. 
Furthermore, GPDs also give information on highly desirable quantities such as the quark 
orbital angular momentum contribution to proton spin [2]. 

Useful constraints on the forms of the nucleon GPDs have been obtained at DESY [10] 
and Jefferson Lab [11]. Typically the GPDs contribute to experimental observables through 
convolutions; they are not directly measurable through those experiments. Thus inputs 
from theories are important. Valuable insights are gained through the impact parameter 
distribution interpretation [7] (see also [8]) and model computations of GPDs [see e.g. [6] 
for a review] . Recently, lattice techniques have first been applied to compute the moments of 
nucleon GPDs directly from QCD [12]. The latest unquenched lattice results are presented in 
[13]. However, due to the limitation of computing power, these calculations often employed 
u and d quark masses heavier than their physical values and employed big momentum 
transfer. Thus, extrapolations in quark mass and momentum transfer are required to obtain 
physical results. In general, the quark mass and momentum transfer dependence have non- 
analytic structures which should be incorporated in the parametrization of the extrapolation 
formulas. Fortunately, chiral perturbation theory (xPT) [14-17], which is an effective field 
theory of QCD, can be applied to extract those non-analytic structures in a systematic, 
model-independent way. 

Recently xPT has first been applied to the computation of hadronic twist-2 matrix ele- 
ments [18, 19], which is related to the moments of PDFs and GPDs through the operator 
product expansions. Many applications have been worked out, e.g., chiral extrapolations 
of lattice data [18-21] including (partially) quenching [22, 23] and finite volume [21] ef- 
fects, GPDs for quark contribution to proton spin [24] , gravitational form factors [25] , pion 
GPDs [26, 27], large Nc relations among nucleon and A-resonance distributions [28] (see 
also earlier work in the large Nq limit [29]), soft pion productions in deeply virtual Compton 
scattering [30-32], SU(3) symmetry breaking in the complete set of twist-2 [33] and twist-3 
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[34] light cone distribution functions, pion-photon transition distributions [35] and exclusive 
semileptonic B decays [36]. The method is also generahzed to the multi-nucleon case [32, 37]. 
There are also earher results derived using the soft pion theorem [38-41] . 

In this paper, we use xPT to study the complete set of nucleon GPDs. We compute the 
leading quark mass and momentum transfer dependence of the moments of nucleon GPDs 
through computing the nucleon off-forward (meaning the initial and final nucleon momenta 
being different) twist-2 matrix elements. Then we apply the xPT results of the moments to 
get insight on the GPDs themselves and their impact parameter space distributions. 



II. NUCLEON GPDS 



The eight nucleon GPDs, H, E, H, E, Ht, Ht, Et, and Et, form the complete set of 
leading twist GPDs for quarks of flavor q [42] . They are defined through off-forward nucleon 
matrix elements of the vector, axial vector and tensor quark bilinears: 



(P'|g(-|nyg(|n) \P) = y dxe— ^[/(P') 



Ht^ + E 



2M 



{P'\Q{~n)i/i^,q[^n) |P) 



U{P): (1) 
U{P) , (2) 



(P'|g(-^n)zn«r;3a"^g(|n) |P) 



{n • P) (r • A) - (r • P) (n • A) 



#(r-A)-/(n-A) , ~ ^ {r -P) - f {n -P) 

-THiJ' h HiX 



2M 



2M 



M2 
U{P) 



(3) 



where M is the nucleon mass, U is the nucleon Dirac spinor with normalization U {P)U{P) = 
2M, pi" = {P + P'Y/2, A^ = P^' - P/", n is a dimensionless light-like vector [n^ = 0], 
^ = —n ■ A/ {2n ■ P) and t — A^. r is a transverse vector satisfying r • n — r • P — 0. We 
have used e°^^^ = — 1 and the notation e^^^*^ — e'^^'^^AaBpC^Ds- Here we have suppressed 
the Wilson lines connecting quark fields locating at — |n and |n to make the nonlocal quark 
operators gauge invariant. These GPDs are functions of x, ^, and t. Both x and ^ have 
support from —1 to -fl. 

The GPDs encode information of ordinary PDFs and nucleon form factors [2]. In the 
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forward limit of A'* —> 0, we have 

H{x,0,0)^h{x), H{x,0,0)^gi{x), Ht{x,0,0) ^ Mx), (4) 

where fi{x), gi{x) and frix) are spin-averaged, hehcity and transvcrsity PDFs, respectively. 
On the other hand, forming the first x moments of the new distributions, one gets the 
following sum rules, 

J dxH{x,tt) = Fi{t) , 
J dxE{x,^,t) = F2{t) , 
j dxH{x,C,t) = GA{t) , 

j dxE{x,^,t) = Gp{t) . (5) 

where Fi and F2 are the Dirac and Pauli form factors and Ga and Gp are the axial- vector 
and pseudo-scalar form factors. There are also tensor form factors associated with the first 
X moments of Ht, Ht and Et, but not Et, because time reversal invariance demands [42] 

J dxET{x,tt) = . (6) 

The most interesting sum rule relevant to the nucleon spin is [2], 

J<i = lJ dxx[H{x, 0) + E{x, 0)] , (7) 

where Jg is the q quark contribution to proton spin in a frame in which the proton has a 
definite helicity. The ^ dependence in the sum rule has dropped out. Since Jg can further be 
decomposed into quark helicity and orbital angular moment contributions, by measuring Jg 
from experiments sensitive to GPDs and measuring the helicity contribution from polarized 
deep inelastic scattering, the quark orbital angular momentum contribution to proton spin 
in principle can be obtained. 

There are also a set of eight gluon GPDs defined as the matrix elements of non-local gluon 
operators. They mix with the isoscalar combination of quark GPDs under renormalization 
scale and transform in the same way as isoscalar quark GPDs under chiral transformation. 
We will first focus on the quark GPDs, and later come back to the gluon GPDs. 
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III. CHIRAL PERTURBATION THEORY 



XPT is a low-energy effective field theory of QCD. xPT makes use of the symmetries 
and scale separation of QCD and allows a model independent description of physics below 
the chiral symmetry breaking scale /^^(^ ^ttF^ ~ 1 GeV), where = 93 MeV is the pion 
decay constant. In this paper, we will compute single nucleon matrix elements associated 
with nucleon GPDs; thus, the relevant scales below /j,^ (light scales) include the pion mass 
TTT-TT ~ 139 MeV and the characteristic momentum p in the problem. The nucleon mass M, 
which is numerically of the same size as /i-^, is treated as a heavy scale as /i^. Thus the 
standard heavy baryon xPT [16] approach is used to systematically disentangle the hght 
and heavy scales. Here the following four small expansion parameters are treated the same 
in the chiral expansion and denoted as 

^ ^ 12^' 12^' M' M' ^ ' 

The physical pion fields (tt^, tt"*", tt") enter the theory through the matrices 

E = e'"/^% u^V^, (9) 

where 

n = . 10 

The relevant terms in the chiral Lagrangian are 

C = ^Tr [ d^E d^E^ ] + A Tr [ m^E^ + h.c. ] + 
8 

+iN^v ■ DN + 2gAN'^ (11) 

where the quark mass matrix rUq — diag{mu,mci) — mj, and we will take the isospin 
symmetry limit m„ = rud — m. The nucleon field N — (p, n)^, v is the nucleon velocity 
and S'^ — ^a'^^^^Vi, is the nucleon spin vector, v ■ S — qa — 1-26 is the axial-vector 
coupling constant. The pion-nucleon couplings arise in Eq. (11) through the vector and 
axial couplings 

= ^ {ud^'v) + v^df^u) , A'' = l- {ud^'v) - v)d^u) , (12) 

and the chiral covariant derivative 

DM ^ (QM + v^) . (13) 



Under a SU{2)l(S) SU{2)r chiral rotation, the hadronic fields in the chiral lagrangian trans- 
form as 

E — > LUR^ , ruq — > LniqR^ , 
N U{x)N , u LvU\x) = U{x)uR^ , 

U{x)A^U^{x) , D^N U{x)D^N , (14) 

such that the chiral lagrangian in Eq. (11) stays invariant under chiral transformation. 
IV. THE VECTOR OPERATORS 



Instead of directly matching the non-local quark bilinear operators to hadronic operators 
in xPT, it is conceptually more straightforward to deal with matching of local operators. 
One can perform operator product expansions (OPEs) to convert the non-local operators 
to the sums of local twist-2 operators then do the matching [18, 19]. The Taylor-series 
expansion of Eq.(l) gives 



(P'|0'"|P) = {n-P)""U{P') 



m+ll 



+ E, 



m+l' 



2M 



(15) 



where 



in-D] q (16) 

is a twist-2 operator dotted by the n^°n'^^ • ■ ■ n^" tensor to project out the symmetric and 
traceless part. The gauge invariant covariant derivative — {D^^ — D^^)/2. Hjn+i(^,t) — 
J dxx"^H{x, ^, t) and Ejn+i{^, t) — J dxx^E [x, ^, t) are the (m -|- l)-th moments in x of the 
GPDs. The nucleon matrix element of has different form factor structures [2, 43, 44]. 
Using the notation of Ref. [44], 



(P'|0"*|P) = U{P') 



rn 



j=0 
even 



—i- 



2M 



[n 



AY {n ■ P)"'-' B^^,,{t) 



m odd- 



U{P) 



(17) 
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The constraints on j and Cm are due to the requirement of time reversal invariance [2] . By 
comparing this equation with Eq.(15), we have 



I m odd ' 



/I m 

even 

/I m 
dxx-E{x,i,t) = Y.{-2iy Bm+iAt) - (-20™+' Cm+mm odd m 

even 

after using the Gordon decomposition. The time reversal invariance demands i^m+i and 
Ejn+i to be even in ^: 

Hm+i{-C,t) = Hm+i{C,t) , Em+i{-C,t) = Em+i{C,t) . (19) 

To apply heavy baryon xPT, we perform the 1/M expansion to Eq.(17) [16]. The leading 
terms are 



{F\0"'\P) = 2N 



m 

J2{(n-Ay{Mn.vr-'^'Em^,,{t) 



j=Q 

even 



where 



+ ze'^"^^ (n ■ A)-'' (Mn ■ vf'-' M^+i,,(t)} 

+ (n.Ar^C„^,(t)L JiV, (20) 



Em+l,j — Am+l,j + ' ^m+l,3 — ^m+lj + -^m+lj ■ (21) 



For convenience, we will work in the Breit frame where = (1, 0) = P^^ /M + 0{1/M), ~ 
(0, a/2) and = (0, A ). The normalization of the Pauh spinor N is iViV = 1 + 0(1/M). 

A. Pionic Vector Operators 

In a similar manner, the pion vector GPD is defined as 

= J dye-'y^^-^-H^{y, e., t)n ■ P^tr [rW] , (22) 

where the isospin operator r" = (1, V), and from now on q is an isospin multiplet [but 
note that the isoscalar quark also contains the s quark contribution]. There is no GPD 



because the pion is spinless. The Taylor-series expansion of the above equation gives 

{AK)m^KP.)) = t) {n ■ P.)™+' tr [tW] , (23) 

where = (P^ + P^) /2 and 

C>^ = gr"^(n-zB)™g, (24) 

and H^,m+i{^7v: = / dyv'^H^iV: Ctt, *)■ 

The pion GPDs are strongly constrained by charge conjugation (C) and isospin symmetry. 
Under C, 

CQrn^-l = ^ (25) 

for a = and 3 [33]. For m even (m = 2k), the above equation implies 

<7r°(P;)|Oo^V°(^0> = , <vr°(P:)|Of k°(^.)> = , (26) 

because 7r° is C even. Furthermore, by isospin symmetry, (7r='=(P^)|OQ*^|7r='=(P^)) = 0. This 
implies 

<2fe+i(^.,i) = 0, (27) 

or, equivalently, 

i/°(y,e.,t) = i/°(-y,e.,t) . (28) 

For m odd (m = 2^ — 1), 

(7r^(P;)|Oo^^-V^-(PO) 
= 2%||^(n-A)^^(n-P.)^'-^'^-^,,(t) + (n-A)^^C-°(t)| , (29) 

where P^ - P^ = A and = -n • A/ (2n • P^) . 

On the other hand, for the isovector case, under C 

= -(7r+(P;)|OrV-'ra) , (30) 

where the equals sign is due to (Og*^"^) oc (T3). This, together with (7r°|03*^~^|7r°) = 0, 
implies 

^2.(^^,0 = 0. (31) 



Thus 



Again, time reversal invariance requires H^{y,^Tj.,t) to be even in ^^r- 

H:{y,C.,t)^H:{y,-C.,t) , (33) 
as shown in Eqs. (29) and (32). 

In matching to the hadronic operators in xPT, it is useful to write 

= + C)^^ , (34) 

where 

0':,R = QR rliin ■ iDTqn , (35) 

and similarly for O^^. qL,R — [(1 T 75)/2]? is the left(right)-handed quark field. The 
distinction between t£ and is only for bookkeeping purposes. We will set t£ = = r" at 
the end. Under a global chiral SU{2)l x SU{2)ji transformation, — > iig^^ and qL Lq^. 
If we demand 

Tl ^ Lr^Lt, ^ i?r^i?t, (36) 

then O'^^^j^^^ will be invariant under chiral transformation. Furthermore, under charge con- 
jugation, E ^ E^, ruq —>■ rrig, and if we demand 

rl - rf , r» ^ r'f , (37) 

then Eq.(25) can be satisfied for any a. 

Using the symmetries mentioned above, we now match to the most general combina- 
tion of hadronic operators with the same symmetries, 

+ + , (38) 

where denotes hadronic operators made of purely pion fields and O^jy denotes hadronic 
operators with nucleon number equal to one. The ellipse denotes operators which do not 
contribute to our SU(2) calculations in single nucleon sectors, such as operators with nucleon 
number equal to two and above or operators with hyperon fields. 
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For a given m,the leading pionic operators in the chiral expansion are 



7-1 9 

OZ^- ^Y. ^rnl.,j (0) {-in ■ ay tr r^E (in • d) Et 



even 



m—j+1 



. I X m-j+l 



+ 



(39) 



where = {d ^ — d^)/2 and the eUipse denotes higher order operators with more powers 
of derivatives or quark mass matrix, m = 0, 1, 2... There is no restriction on the value of m 
[18, 19]. For a = (the isoscalar case), it is easy to see that the contribution with m even 
vanishes — a consequence of charge conjugation. Using 



-in ■ d) 



2m 



(in - a)^"Et 



-m ■ a) 

{—m ■ a) 



E (in ■ d) E^ 

{in ■ (9S) (in ■ {in ■ dH^) 



2n-2 



the operator in Eq.(39) can be rewritten as 



fc-i 



1=0 



[m 



/ ^\2k-2l-2 

dE) (in ■ d) {in ■ dE^) 



+ •••, (40) 



which is the operator constructed in Ref. [26]. The prefactors in Eq.(40) are chosen such 
that Eq.(29) is reproduced in the leading order in the chiral expansion subject to constraints 
that relate 6*2^° to A^f^i- The constraints come from the (n ■ Pj,) (n ■ P^) = (n ■ P-j^Y (1 — ^^) 
factor in {n'{P;^)\Ol''-^\n^ {P^)) which makes H^{y,^l = ±1,0) = [this was also observed 
in Refs. [38, 39]]. This property, however, does not persist at higher orders. At the next-to- 
leading order (NLO), there several sets of counterterms: 



^0,7r,l — "2fc 

k-1 



{—in • d)^^ tr [Em^ -|- m^E^] , 



k-l 



, 2k— 21— 2 

Oo5rT2 = Yl ^2k,2lt^ [^H + ^<?^T {-^^ ■ {^^ ■ (^^ • ^) (^^ • 



1=0 
k-l 



(/ 2k— 21— 2 

{m ■ aE) (in ■ d) {in ■ ^Et) 



{in ■ 9E) (in-d 



2k-2l-2 



{in ■ 9Et) 



(41) 



It is O^V/ that makes H^{y, ±1, 0) non- vanishing. 
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B. Nucleon Vector Operators 



In matching to nucleon operators, it is convenient to define 



(42) 



such that under a chiral transformation, r",^ — > U{x)T^_^U{xy with T^(i) transforms as in 
Eq.(36). The leading operators contributing to ^m+i 

^vnAS 



+ 2M{Mn-v)"'N 



<?i(e,0)n.. + C^S(e,0)^ 



"(0) 



M 

SnAS 



2M 



T,? iV + 



(43) 



where -Em+i = Hm+i + Em+i- The prefactors are chosen such that Eq.(15) is reproduced in 
the leading order in the chiral expansion and a superscript (0) denotes the chiral limit value. 
The Crn+i and D^+i operators will not contribute at the order we are working (NLO). The 
Hm+i and E'm+i operators in Eq.(43) can be further written as 



O 



a,N 



2N 



m 

5]{(n.A)^(Mn.t;)-^-+^KS.(0) 



+ ie 



j=0 

even 

vnAS 



{n ■ Ay {Mn ■ vr-' mZUo)] 



+ (n-A)-+^CS(0) 



m odd- 



(44) 



where we have replaced the total derivative operator —in • D by n • A. The operators 
with derivatives acting only on t^_,_ or do not contribute to nucleon GPDs by direct 
computation. 



C. Form factor results of vector twist-2 matrix elements 

Now wc present the leading chiral corrections to the form factors of the vector twist-2 
operators defined in Eq.(20). We will insert powers of £ to keep track of the chiral expansion. 
One should set £ = 1 when using these results. The results of £'^^i q(0) reproduce those 
of the forward twist-2 matrix elements in Refs. [18, 19, 28], while -E'J'o(^) and M^o(i) agree 
with the electric and magnetic form factor results in xPT [17]. 
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(a) (b) (c) 




FIG. 1: Diagrams contributing to leading chiral corrections of nucleon off- forward twist-2 matrix 
elements. The squares denote insertions of hadronic operators. The solid lines are nucleon fields 
and the dashed lines are pion fields. The diagram (f) denotes the wave function renormalization 
contribution. The leading order pion pole diagram is shown in Fig. 2. 

For the isoscalar {a = 0) form factors^, 



= (l + ^n)MS,2.(0) + ^W+iI^2i\2.W 
+^'7^+i,2fc< + ^'C+i,2fci + Oie') , 
C^^iit) = 0(0) + e6r^,,,^,Vg_,, it) + 0{e') , (45) 

where = 0, 1, 2 . . . The C contribution is from the wave function renormalization (Fig. 1(f)) 
and the loop diagrams with one insertion of the nucleon operators in Eq. (44) (Figs. 1(b) 
and (c)). The V contributions are from the loop diagrams with one insertion of the pionic 
operators in Eq. (44) (Figs. 1 (d) and (e)). 

For the -E^+i 2fc(^) form factor, the one-loop contributions from Figs. 1(b) and (f) cancel 
each other while the diagrams in Figs. 1(d) and (e) are of higher order and Fig. 1(c) 



^ There is no restriction on the range of m in this ChPT calculation. We are only interested in the to^ and 
t dependence of the form factors, thus we keep track of the expansion in m.^ and t using the parameter e 
and count x = 0{e'^). There is no intrinsic difference between operators of different m — they all transform 
in the same way under a chiral rotation [18, 19]. 
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vanishes. At 0(£^), -E^+i,2ifc ^-Iso receives analytic contributions proportional to or t 
from counterterms (Fig. 1(a)) except for the charge operator £'°q(0). a^+i^2k ^^'^ Pm+i,2k 
are independent of m^, t and the renormalization scale /i. 
For the gfe form factor, 

and 

,2 *= /•! /I \2; 



1^' 



^2^+2,2. (^) = E ^2^+2,2.-2. (0) dy {I + l) (2Z + 1) " I/) ^(y)^ log . 

(47) 

where the integration variable y arises from the Feynman parametrization and m{y) = 
^/ml — y{l — y)A'^. The counterterm contributions 7^+1 2fe a-nd Cm+i,2fe depend on but 
not on uIt^ and t. The /i dependence of j!^+i^2k Cm+i,2fe cancels the /i dependence from 

As for the C^+i form factor, it receives non-analj^ic contributions from the Fig. 1(d) 
diagram at 0{e) and no contribution from analytic counterterms until 0{e'^). 



-^^2 it) = ^ |e ^2^2,2.-2. (0) I' dy(^l-y 

+ ^2^2(0) f'dy 
Jo 



21+2 



2 



m(y) 



^ - 4 (/ + 2) m(y) 



4m(y) 



(48) 



_m(?/) 

By setting {m,k) = (1,0) and using ^2,0(0) — — 46*2 '° (0) = (^)7r' results in Eq.(45) 
reproduce those of Refs. [24] and [25] . 

The leading chiral corrections for the isovector {a — 3) form factors are 



-^m+l,2fe(^) - (1 + ^^(^e) -E^m+l,2fc(0) + ^^<^m,2fe'P2fc+l,2fe (0 

(1 - (5™,o(^fc,o) a^+i,2fc"^' + £''Pi+i,2ki + C(^^) , 
M^+,,2kit) = M'J:%,{Q)+e5rn,2kVll'l,^,,{t) + 0{e') , 

^^^-i(i) = (1 + e'Cl) CS(0) + s^^^+.m^ + e^l^.t + 0{e') . (49) 
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Here 



and 



r)3,M 



— 



— 



ml 



(Sgl + l)ml 

{Sgl + l)ml ml 
(47rF.)2 



(50) 



- V A^'^ ((]) ^^^^ 



=0 " 

" ^) " + (5 + 4/) m(y)^] + m{y f] log ^ , 



W = - E ^2^+1,2.-2. (O)^J^'ciy I (1 + 20 Q 

Here (0) is the number of u quark minus the number of d quark in a tt"*" meson, ^4^0 (0) = 

operator -E'f_o(0) renormahzed. Unhke the isoscalar case, 

' . V 0{e^), respectively. 



= 2; thus the charge wj^v 
the chiral corrections to M^^^ 2fe C^+i start at 0(£) 



V. THE AXIAL OPERATORS 



The Taylor-series expansion of Eq.(2) gives the form factors of off-forward nucleon axial 
twist-2 matrix elements 

m 

{P'\Ota\P) = E^(^')f#75(n-A)^(n-P)'"-^iV,(^) 
1 



j=0 

even 



+ 75^(n-Ar^(n-P)'"-^SV,.(t) 



t"C/(P) , 



where 



For m — 0, it reduces to the nucleon axial current matrix element 



(52) 



(53) 



(54) 
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with ^1^0 = Ga and Bi^ — Gp. The form factors are related to the moments of axial GPDs 



as 



/I ^ m 

even 

/I m 

even 



(55) 



Similar to the vector twist-2 operators, O™^ is matched to pionic and nucleon operators 

OT,a - + 05™.,iv + • • • ■ (56) 

Since O^^, = O™^ — O™^, the matching to leading pionic operators is similar to that of 
Eq.(39) with a different relative sign, 



j=0 



, m-i+1 



+r^E^ fin • a j E 



+ ••• (57) 



Parity conservation governs that the axial operators match onto pionic operators with odd 
number of pions; thus diagrams in Figs. 1(d) and (e) could not contribute. Instead, the 
pion pole diagram in Fig. 2 gives the leading contribution from insertion of O™^^. In this 
diagram, a necessary input is the tt — > matrix element of O^^^ .^^ which is related to the 
(m + l)-th moment of pion light cone distribution function 

(7r''(A)|0- |0) = tS^,F^ i-n ■ A)™+^ {z^^ /2- . (58) 

where z — 1 — 2x', 

{z^)^= f dx'{l-2xT<t>A^') , (59) 
and where 0,^ {x') is the pion light cone distribution function with the normalization {z^)^ = 
dx'(t)T^ {x') = 1. {z"^)^ vanishes for odd m due to charge conjugation. In Ref. [33], it was 
shown that the leading non-analytic quark mass corrections to {z'^) can all be absorbed into 
F^. Thus {z^)^ is purely analytic at 0{e^) [33] but not analytic at higher orders [34]. Using 
the operators constructed in Eqs. (39) and (57), our leading-order (LO) result reproduces 
that of Ref. [38] derived from the soft pion theorem, 

Hr\x,i,^l,t^Q)^\4>.(^-^^ ■ (60) 
15 



FIG. 2: The leading order pion pole diagram contributing to the nucleon pseudo-scalar form factors. 



In heavy baryon xPT, the leading axial nucleon operators in the matching are 

m 

OZ,N = E f2n .S{n. Ay {Mn ■ t;)-^' ^:?1,(0) 

CiV + ■ ■ ■ , (61) 



j=0 
even 



where 

Em+l,j = ^m+l,j 5 ^m+l,j = , '"^ = + B^+IJ ■ (62) 

For the isoscalar case, the leading chiral corrections include diagrams of Figs. 1(a), (b), 
(f): 

E^^i,2kit) = (1 + e'Cl,) eT^,^,,{0) + e'al^,^,,ml + 5^^+1,2.^ + 0{e') , (63) 

The contribution of E^'^^ ^m+i 2k from 1/M^ corrections of Fig. 1(b) type diagrams 
[45]. 

For the isovector case, the leading non-analytic contribution of is from Figs. 1(b), (c) 
and (f). 

Ei+,,2kit) = (1 + e'Clj) ^S,2fc(0) + ^'5^+1,2^^^' + ^'£+1,2.^ + Oie') , (64) 

where 

{2gl + l)ml ml 
" (4vrF.)2 • ^^^) 

The /i dependence in Eqs. (63) and (64) is absorbed by the counterterms a'^^^ 7m+i 2k 
and 5^+i,2fc- 

In contrary, the leading contribution to is non-analytic [(9(£:~^)] arising from the pion 
pole diagram shown in Fig. 2. At one loop, we obtain 

M^+1,2. W = 5^,2kGp{t) {z'')^ /2'' + mS,2.(0) + Oie) , (66) 
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where {z^^)^ is analytic at 0{e'^) as mentioned above and the pseudo-scalar form factor [46], 



Gp{t) 



where dis is a counterterm and (r^) is the square of nucleon axial charge radius. 



(67) 



VI. THE TENSOR OPERATOR 



The tensor operator 



has the nucleon matrix element [44] 



3=0 

{n-P) (r-A)-(r-P) (n • A) 

m+l,j 



J even 



+ 



J even 



i{r-A)-f{n-A) T,a ... 



2M 



J even 



+ 



{r-P)-f{n-P) ~ 



M 



j odd 



t'^U{P) . 



The form factors are related to the moments of nucleon tensor GPDs as 

HT,m+l 



E' 



T,m+1 



T,m+1 



ET,m+l 

After the 1/M expansions, Eq. (69) becomes 



/I m 
j=0,even 

/I m 
dxx^^Erix^U)- E (-20^- S^+i,,.(i) , 
j=0,ei)en 

/I m ^ 

dxx-^T(a;,e,t)= 5^ (-20^- i^+,,,.(t) , 
j=0,even 

/I _ _ 

■1 n 



j=0,odd 



(P'lOJpjP) = J2 '^MN{P') {n ■ AY {Mn ■ v) 

3=0 



m-j 



2i6-^M^fi,.-(0 



J even 



+ 



{n-v){r-A) T,a 



M 



-^m'+lj(^) 



^vrAS 



J even 



i{n-v)e T,a 



M 



(n • A) e*"^^^ - (r ■ A) e'^"^^ ta , . 
+ i- ^2 <'tiAi) 



J even 



t"N{P) , 



(68) 



(69) 



(70) 



(71) 
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with 

Under charge conjugation, 

CO??,,C-^ = (-1)"^+' 0?^a (« = 0, 3) . (73) 

By similar arguments to the vector operators, the moments of pionic tensor GPD can be 

parametrized as 

{n'{Pi)\q{in^r,an{in ■ Df^+'q^iP^)) 

- 25., ^ [(n • (r • A) - (n • A) (r • P^)] {n ■ Af {n ■ P^f-''^' El%^,,{m) 
1=0 

{7r\P:)\qT'{tn,r,a''ni^n ■ Df^qlw^ (P^)) 

= 2i6''^' E ^ [(^ • ^0 • A) - (n • A) (r • P^)] {n ■ Af {n ■ P^''"'' E^ifiaii^) (75) 
1=1 

where we have inserted a heavy scale M which is of order fiy, to make E^j dimensionless. 

The matching procedure is similar to the vector and axial vector cases which we will not 
repeat here. We just point out one main difference. O'^^ can be decomposed as 

C^a = ?LT^i^^%r,a'^^(^r^ • D)'"g« + qRT^Ltn,r,a''''{tn ■ D^qL . (76) 

We will set = r^^ = at the end. Oj?^ will be invariant under a chiral rotation if we 
demand 

riR^Lr^RR^, t^l^Rt^lL^. (77) 

Thus instead of using r^_|_, we use — \ {u^rf^^u^ ± UTp^iu) which transforms as — > 
U{x)t'^U{x)^ under a chiral rotation to construct the nucleon operators [19]. 
For the isoscalar form factors, 

<ii,2.(^) = (l + ^n/)M:2°L(0)+^^5^'^,2X + ^'ei,^^^ 

EZiAt) = ^SSL(o)+s^e,2X+-^£:M.^ + 0(^^) , 

^m+l,2A:+l(^) = (^ + ^^^m) ^nf+l^k+li^) + ^^'^m,2fe+l'^fe+2,2A;+l (^) 

<i,2.(^) = (1 + e'Cl,) W^f^l,{0) + e^^C°,M:S,(0) + e^S^,,^aZ2,.it) 

+e'^Zi,2kml + e'~XZi,2,t + 0{e') , (78) 

18 



where 

Similar to Eq. (63), the contribution of M'^^^\^, to W[^^^\j^ is from corrections of Fig. 

1(b) type diagrams. 

Similarly, for the isovector form factors. 



<i,2.(^) 

-^m'+l,2fc(^) 
^m+l,2fc+l(^) 



(1 + s'C'^) m2(°),(0) + e'all,^,,ml + £^^1,2.^ + (^(^') , 



zii,2dt) = w^^^i:^.(o)+s5„.,2.7^t:i,2,(^)+(^(e^), 



,T,3(0) 



3,W 



(80) 



where 



and 



_ 2gl + 1/2 



2 1 

— — — In — 



,2 ' 



(81) 



^2fcfl,2feW = - 32^2^2 ^2k+l2k-2li'^) ^ 0^?/ Q " Z/^ 

X {^1 [2m2 - (5 + 4Z) m(y)'] - m(?/)'} log ^ 

A* 



(82) 



It is interesting that the leading chiral corrections of the tensor and vector matrix elements 
are similar. They obey the relations. 



r°'^ (t) - 

■^2k+2,2k+l\^) — 



^fe+l,2fc(^) 



O -^2^+2, 2A: 



"^2^+1,2*; (^) 



-'2fc+l,2fel''J ~ ° -^2fe+l,2fe I'-J 



(•^'+-'-)^2fe+l,2fe-2i^-^2fc+'l,2fc-2I 



,2fc-2i~*-^2fc+'l,2fc-2i 



(83) 
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VII. GLUON GPDS 



The gluon GPDs are defined as [42, 47] 



2M 



2M 



n,n,r^MP'\ F^^ \P) 
= J dxe-^^""--^ [(n • P) (r • A) - (n • A) (r • P)] 

(n • P) (r • A) - (r ■ P) (n ■ A) 



xU{P') 



M2 

{r-T)-f{n-T) 
2M 



(84) 



where the Wilson lines are also suppressed. The Taylor-series expansion in z of the above 
equations gives rise to relations between the gluon twist-2 matrix elements and the moments 
of gluon GPDs: 

Fg,m+i{i,t) = j dxx^Fg{x,U) , (85) 

where Fg{x,^,t) denotes a generic gluon GPD. Here, unlike the quark case, m — 
1,2, . . .without m = 0. This is because the right-hand sides of Eq. (84) are of the form 
xFg instead of Fg. Therefore Fg can only be determined up to a function of the form 
X{^,t)S{x) using the above definitions. However, the first moments in x of gluon GPDs can 
be probed by non-local gauge invariant operators as introduced in Ref. [48]. 

The local gluon twist-2 operators and isoscalar quark twist-2 operators all transform in 
the same way — as singlets — under chiral transformation. The gluon twist-2 operators match 
onto the same set of hadronic operators as the isoscalar quark twist-2 operators with different 
prefactors. Thus our previous results for the moments of isoscalar quark GPDs can be easily 
converted to moments of gluon GPDs. 
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VIII. xPT CONSTRAINTS ON THE NUCLEON GPDS 



In this section we use the moments of nucleon GPDs calculated above to shed light on 
the GPDs themselves. We first pay attention to the non-analytic 0{e) corrections before 
making transit to the impact parameter distributions of GPDs. There are no counterterms 
at 0{e), therefore those corrections are clean predictions of xPT. Here are several interesting 
cases of this type: 

1) At 0{e), and only receive chiral corrections from the "D terms" which only 
depend on x/^ and t [39], 



5E'{x, e, t) = -SH\x, e, t) = t)e (^1 
with D{z,t) = -D{-z,t). From Eqs.(18), (45) and (48), 



(86) 



SH'^^.iU) = j ^ dxx^5H\x,i,t) = {-2ir^'Cm+,{t) , (87) 

where Cm+i{t) = for m even. Then one can prove 5E^ and 5H^ have the functional form 
in Eq.(86). 

2) E^, E^ and receive 0{e) contributions which are predictable. We will discuss their 
impact parameter distributions in the next subsection. 

3) One can show that Eq.(66) imphes 

E'ix, e, t) = -^^Gp{t) + 0{e') . (88) 
This result coincides with those derived in Refs. [40, 41] after using 0^(a;') = 0^(1 — x'). 



A. Impact Parameter Distributions 

In the limit of ,^ — (and t — A^^), the Aj^ Fourier transformation of the GPDs 
has the interpretation of simultaneous measurement of the longitudinal momentum and 
transverse position (impact parameter) of partons in the infinite momentum frame. The 
impact parameter dependent parton distribution for a generic nucleon GPD F is 

J^{x, b^) = / ^^e^^^-^-F{x, 0, -Ai) . (89) 
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We are interested in the small t {t <^ /i^), or large b±{b » 1/a*x) region where xPT is 
reliable, so we expand F in t, 



F{x, 0, t = -Al) = F{x, 0, 0) - dtF{x, 0, 0) + ■ ■ ■ , (90) 

which leads to 

J^{x, h^) = F{x, 0, 0)S^ (bx) + dtF{x, 0, 0) Vl^' (bx) + • • • . (91) 
The F term gives the b± integrated distribution 

J (Ph±T{x, bj_) = F{x, 0, 0) , (92) 

while the dtF term is related to the averaged 6^ of F as a function of x, 

_ J d'h^blTjxM _ d,F{x,0,0) .p.. 
J d''h^J^(x,bj_) F(x,0,0) ■ ^ ^ 

The terms with higher derivatives on the delta function give higher moments of b^. In 
comparison, the charge radii of the electroweak form factors defined in Eq. (5) are 

_ JdxdtFix,0,0) 
\ JdxF(x,0,0) ■ ^ ^ 

They also constrain the functional form of F. 

In the following paragraphs, we will extract from the model-independent results of twist-2 

matrix elements we obtained above to obtain F{x, 0, 0) and {b\)-p- The GPD E^{x, ^, t) has 
a special ^ ^ limit: 



E^{x,0,t) = 6{x)Gp{t) +0{e^) 
ml —te 



where wc have used Eqs.(88) and (67), {x') = for x' < or x' > 1, and dx'cj).,^ {x') = 1. 
This yields {b\)^z = at x = 0. 

For the other GPDs, it is convenient to express F{x, 0, 0) and dtF{x, 0, 0) as (except 
which we will discuss later) 

F(x,0,0) = a^(x)+£6^(a;)^ + £V(x,//)^ + £2d^(x)^logf^) +C»(£3) , 

(96) 



dtF{x,Q,Q) = 4- 



-e^ix)^ + Mx, /.) + gp{x) log ( ^ ) + 0{e) 



X 

.2 
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where all the pref actors are m.^ independent and the // dependence of and fp cancel the 
dependence in the logarithms of the dp and gp terms, respectively. 
We will show that those prefactors could have 5 function structures. Note that instead 
of using 5 functions with zero widths, one can use regularized delta functions with finite 
widths as well. For prefactors with odd powers of e (defined in Eq.(8)), the widths of the 
delta functions should be ~ £ ~ m7r///x or smaller. For prefactors with even powers of £, 
the widths of the delta functions should be ~ £^ ~ (m^///;^)^ or smaller. To demonstrate 
this, we use the regularized delta function — (Ax)"^ 9{/S.x — x)9{x). The x'^ moment 

of e^5/^x{x) gives 

{x"") = / dxe''SAx{x)x"' = ^ ; . (97) 



/o m+1 
When m > 0, the change of {x"^) due to regularization is 0{e'^ (Ax)"*). Thus the change to 

the parton distributions due to regularization is 0{e'^Ax). These contributions should be ab- 
sorbed by matrix elements of higher order operators in ChPT. However, in ChPT, the higher 
order operators are all analytic in light quark mass rriq (oc m^). Thus, the regularization is 
sensible only when (Ax)"* is e to an even power. 

1) For the isovector GPDs E^, Ej, and Hj,, their 0{e) contributions are non- vanishing. 
From Eqs.(18), (21) and (49), the x^ moment of E^{x, 0,t) vanishes except for m = : 

5El,^,{0,t) ^ e5m,oVl^ (t) . (98) 

This implies 

bE3{x) =bE3S (x) , eE3{x) = eE35 (x) . (99) 
where 6^3 = — 126^3 = —g\M/FT^ — —16.0. Thus we have 

d^h^8^{x, b±) - aE3{x) + ebp^S {x) — + 0{e^) , 

^ e ji^m^aE^ix) ^ ^ ^ ^ 

Note that j dxaE^i^x) is the isovector nucleon anomalous magnetic moment in the chiral 
limit. The charge radius square of the isovector nucleon anomalous magnetic moment is 

e ji^rriT^ j dxttE^ix) ^ ^ ' i^^^) 

The behavior of E^{x, 0, t) is similar to that of E^{x, 0, t) with bE3{x) oc 6 (x), 6^3(2;) oc 
5{x), and so on. Also, H^{x,0,t) = -E^{x,0,t)/2 at 0{e) from Eqs. (70), (72) and (80). 

23 



2) For the other GPDs, the 0{e) contribution to F{x,0,t) vanishes such that bpix) — 
epix) — 0. The non-analytic logarithmic terms dpix) and grix) can be predicted in xPT. 
If grix) is non-vanishing at certain x, then the corresponding (fe^)jr diverges in the chiral 
limit. Here we list the results based on the moment computations. 

a) Spin-averaged quark GPDs: The first moment of H{x, 0, 0) gives the total number 
of quarks of certain flavor in the nucleon which is independent of the quark mass. Thus 
J dxcnix) — J (ixdnix) — 0. For the isoscalar combination there are no chiral loga- 
rithms at 0{£^), 

dffoix) ^ gHo{x) = . (102) 

For the isovector combination H^, there are chiral logarithms at 0{e'^) except for the flrst 
moment in x. This implies 



dH^{x) = -{^g\ + l) 



aHi{x) — 5{x) I dxaH^{x 



gHs{x) = . (103) 



The isovector GPD has been discussed in Eq.(lOO). The isoscalar combination E^+H°, 
which satisfies the sum rule in Eq.(7), yields 

dEo{x)+dHo{x) = -3g\ [aEo{x) + aHo{x) + {x)^S' (x)] , gEo{x)+gHo{x) oc S' (x) , (104) 

where S' {x) — dS {x) /dx. This result is derived from Eq.(45) and we have dHo{x) — from 
Eq.(102). 

b) Quark helicity GPDs: is not defined at ^ = 0. For 



And for H, 



and 



d^oix) = -3g'^^aEo{x) - g^ajjoix) , gEo{x) = . (105) 



djjo{x) = -3g\ajjo{x) , gfjoix) = . (106) 



dmi^) = - (2^A + 1) «i?3(a;) , g^six) = . (107) 

c) Quark transversity GPDs: Ex = when ^ = 0. and E^{x, 0, t) are already discussed 
following Eq.(lOl) above. 

dRoi^) = -^9\aH°ix) , 9Hoix) = , (108) 
dndx) = - ( 2c/i + ^ ) ttH^ix) , gndx) = , (109) 



2^ 
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aEO^{x) + ^aHo{x) + 2E^;,^'\0)S'{x) 



gEo{x) oc S'{x) , 



(110) 
(111) 



-^2,b"' (0)5' (a:) 

3) In Ref. [49], Strikman and Weiss argued that the (fei)-^^ of the gluon GPD Hg{x, b±) is 
proportional to at small x {x < m.^/M) and is proportional to 1/M^ otherwise. While 

this behavior is qualitatively similar to E^{x,b±) [also E^{x,b±) and Hj.{x,bj_)] described 
in Eq.(lOO), for Hg{x, b±) the result is still inconclusive in our calculation. We confirm that 

(bl) 



0{1/M'^) or (9(1//^^) when x ^ from the result for the isoscalar moments in 
Eq.(45) but with m > 1. However, the possible S{x) contribution cannot be probed by the 
matrix element defined in Eq. (84). We will defer the investigation to a later publication. 

Finally, Eqs.(104), (110) and (111) imply that x {VPg + f^), xE^g and xE^g have {b\) oc 
5{x) log (m^) which diverges in the chiral limit. 



IX. CONCLUSIONS 



Using heavy baryon xPT we have studied the leading chiral corrections to the complete 
set of nucleon GPDs. We have computed the leading quark mass and momentum transfer 
dependence of the moments of nucleon GPDs through the nucleon off-forward twist-2 matrix 
elements. We have also applied these results to get insight on the GPDs and their impact 
parameter space distributions. 
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